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1 Stellar Explosions at the Edge of the Visible Universe 

Gamma-Ray Bursts (GRBs) are believed to originate in compact remnants 
(black holes or neutron stars) of massive stars. Their high luminosities make 
them detectable out to the edge of the visible universe |231 IT2j . GRBs offer 
the opportunity to detect the most distant (and hence earliest 1 ) population 
of massive stars, one star at a time. In the hierarchical assembly process of 
halos which are dominated by cold dark matter, the first galaxies should have 
had lower masses (and lower stellar luminosities) than their low-redshift coun- 
terparts. Consequently, the characteristic luminosity of galaxies or quasars is 
expected to decline with increasing redshift. GRB afterglows, which already 
produce a peak flux comparable to that of quasars or starburst galaxies at 
z ~ 1 — 2, are therefore expected to outshine any competing source at the 
highest redshifts, when the first dwarf galaxies have formed in the universe. 

Preliminary polarization data from the Wilkinson Microwave Background 
Probe ( WMAP) indicates an optical depth to electron scattering of ~ 17 ±4% 
after cosmological recombination This implies that the first stars must 
have formed at a redshift z ~ 20 |33l 1111 1131 EHj and reionized a substantial 
fraction of the intergalactic hydrogen around that time. Early reionization can 
be achieved with plausible star formation parameters in the standard ylCDM 
cosmology; in fact, the required optical depth can be achieved in a variety 
of very different ionization histories (since WMAP places only an integral 
constraint on these histories ^H]). One would like to probe the full history of 
reionization in order to disentangle the properties and formation history of 
the stars that are responsible for it. GRB afterglows offer the opportunity to 
detect stars as well as to probe the ionization state |Kj and metal enrichment 
level of the intervening intergalactic medium (IGM). 

1 Observational cosmology resembles archeology. By probing deeper into the uni- 
verse, one reveals layers of it that are more ancient (due to the finite speed of 
light). 
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GRBs, the electromagnetically-brightest explosions in the universe, should 
be detectable out to redshifts z > 10 (23] E|- High-redshift GRBs can be 
easily identified through infrared photometry, based on the Lya break induced 
by absorption of their spectrum at wavelengths below 1.216 /im [(1 + z)/10]. 
Follow-up spectroscopy of high-redshift candidates can then be performed on 
a 10-meter-class telescope. There are four main advantages of GRBs relative 
to traditional cosmic sources such as quasars: 

• The GRB afterglow flux at a given observed time lag after the 7-ray trigger 
is not expected to fade significantly with increasing redshift, since higher 
redshifts translate to earlier times in the source frame, during which the 
afterglow is intrinsically brighter |12j . For standard afterglow lightcurves 
and spectra, the increase in the luminosity distance with redshift is com- 
pensated by this cosmological time- stretching effect. 

• As already mentioned, in the standard ylCDM cosmology, galaxies form 
hierarchically, starting from small masses and increasing their average mass 
with cosmic time. Hence, the characteristic mass of quasar black holes 
and the total stellar mass of a galaxy were smaller at higher redshifts, 
making these sources intrinsically fainter |32| - However, GRBs are believed 
to originate from a stellar mass progenitor and so the intrinsic luminosity 
of their engine should not depend on the mass of their host galaxy. GRB 
afterglows are therefore expected to outshine their host galaxies by a factor 
that gets larger with increasing redshift. 

• Since the progenitors of GRBs are believed to be stellar, they likely origi- 
nate in the most common star-forming galaxies at a given redshift rather 
than in the most massive host galaxies, as is the case for bright quasars 
5 . Low-mass host galaxies induce only a weak ionization effect on the sur- 
rounding IGM and do not greatly perturb the Hubble flow around them. 
Hence, the Lya damping wing should be closer to the idealized unper- 
turbed IGM case and its detailed spectral shape should be easier to 
interpret. Note also that unlike the case of a quasar, a GRB afterglow can 
itself ionize at most ~4x 10 4 E 5 iMq of hydrogen if its UV energy is E 5i 
in units of 10 51 ergs (based on the available number of ionizing photons), 
and so it should have a negligible cosmic effect on the surrounding IGM. 

• GRB afterglows have smooth (broken power-law) continuum spectra un- 
like quasars which show strong spectral features (such as broad emission 
lines or the so-called "blue bump" ) that complicate the extraction of IGM 
absorption features. In particular, the continuum extrapolation into the 
Lya damping wing (the Gunn-Peterson JH] absorption trough) during 
the epoch of reionization is much more straightforward for the smooth UV 
spectra of GRB afterglows than for quasars with an underlying broad Lya 
emission line ^ |S] . 

Although the nature of the central engine that powers the relativistic jets 
of GRBs is unknown, recent evidence indicates that GRBs trace the formation 
of massive stars [221 CUB Ell IE] j and in particular that long-duration GRBs 
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are associated with Type Ib,c supernovae |2H]- Since the first stars in the 
universe are predicted to be predominantly massive [TJ |Sj , their death might 
give rise to large numbers of GRBs at high redshifts. 

2 GRB Afterglows at High Redshifts 

2.1 Motivation 

In difference from quasars of comparable brightness, GRB afterglows are 
short-lived and release ~ 10 orders of magnitude less energy into the surround- 
ing IGM. Beyond the scale of their host galaxy, they have a negligible effect 
on their cosmological environment 2 . Consequently, they make ideal probes of 
the IGM during the reionization epoch. Their rest-frame UV spectra can be 
used to probe the ionization state of the IGM through the spectral shape 
of the Gunn-Peterson |18j (Lya) absorption trough, or its metal enrichment 
history through the intersection of enriched bubbles of supernova ejecta from 
early galaxies^HI- Afterglows that are unusually bright (> lOmJy) at radio 
frequencies should also show a detectable forest of 21 cm absorption lines due 
to enhanced HI column densities in sheets, filaments, and collapsed mini-halos 
within the IGM HH1 

Another advantage of GRB afterglows is that once they fade away, one 
may search for their host galaxies. Hence, GRBs may serve as signposts of 
the earliest dwarf galaxies that are otherwise too faint or rare on their own 
for a dedicated search to find them. Detection of metal absorption lines from 
the host galaxy in the afterglow spectrum, offers an unusual opportunity to 
study the physical conditions (temperature, metallicity, ionization state, and 
kinematics) in the interstellar medium of these high-redshift galaxies. 

2.2 Afterglow Flux at the Lye* Wavelength 

Figure 1 shows the expected spectral flux from GRB afterglows as a function 
of observed time after the GRB trigger for a sequence of redshifts (assuming 
typical model parameters |o])- The flux is calculated at the rest-frame Lya 
wavelength, where intergalactic HI should produce the Gunn-Peterson [T%] 
trough prior to reionization. As already mentioned, the flux does not decline 
dramatically with increasing redshift since cosmic time stretching counteracts 
the luminosity distance increase. 

The expected spectral shape of the Gunn-Peterson trough prior to the 
final overlap phase of HII regions in the IGM, is shown in Figure 2. If 
GRBs trace the typical sites of star formation, then most of them should be 

2 Note, however, that feedback from a single GRB or supernova on the gas confined 
within early dwarf galaxies could be dramatic, since the binding energy of most 
galaxies at z > 10 is lower than 10 51 ergs 
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Fig. 1. GRB afterglow flux as a function of time since the 7-ray trigger in the 
observer frame (taken from The flux (solid curves) is calculated at the redshifted 
Lya wavelength. The dotted curves show the planned detection threshold for the 
James Webb Space Telescope (JWST), assuming a spectral resolution R = 5000 with 
the near infrared spectrometer, a signal to noise ratio of 5 per spectral resolution 
element, and an exposure time equal to 20% of the time since the GRB explosion 
(see |http:/ /www.ngst .stsci.edu/nms/main~ \ . Each set of curves shows a sequence 
of redshifts, namely z = 5, 7, 9, 11, 13, and 15, respectively, from top to bottom. 

hosted by dwarf galaxies at high redshifts. Low-mass hosts would perturb only 
weakly the surrounding IGM (radiatively through their ionizing flux, grav- 
itationally [through cosmological infall, or hydrodynamically through their 
winds), in difference from the massive host galaxies of the brightest quasars 
at the same epoch 0] [5] . Consequently, the spectral shape of the Lya trough 
(which is proportional to e~ T ( x \ where r(A) is the optical depth as a func- 
tion of wavelength) is simpler to interpret for GRB afterglows than it is for 
quasars. The optical depth as a function of wavelength in the Lya damping 
wing is linearly proprtional to the mean neutral fraction of the IGM, xhi, 
since its amplitude is normalized by the optical depth at the Lya resonance, 
~ 6.75 x 10 5 xm[(l + z)/10] 3 / 2 [230. 
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Fig. 2. Expected spectral shape of the Lya absorption trough due to intergalactic 
absorption in GRB afterglows (taken from pj). The spectrum is presented in terms 
of the flux density F v versus relative observed wavelength AX, for a source redshift 
z = 7 (assumed to be prior to the final reionization phase) and the typical halo mass 
M = 4 x 10 8 Mq expected for GRB host galaxies that cool via atomic transitions. 
Top panel: Two examples for the predicted spectrum including fGM HI absorption 
(both resonant and damping wing), for host galaxies with (i) an age ts = 10 7 yr, 
a UV escape fraction / csc = 10% and a Scalo initial mass function (IMF) in solid 
curves, or (ii) ts = 10 8 yr, / csc = 90% and massive (> WOMq) Pop III stars in 
dashed curves. The observed time after the 7-ray trigger is one hour, one day, and 
ten days, from top to bottom, respectively. Bottom panel: Predicted spectra one 
day after a GRB for a host galaxy with ts = 10 7 yr, / csc = 10% and a Scalo IMF. 
Shown is the unabsorbed GRB afterglow (dot-short dashed curve), the afterglow 
with resonant IGM absorption only (dot-long dashed curve), and the afterglow with 
full (resonant and damping wing) IGM absorption (solid curve). Also shown, with 
1.7 magnitudes of extinction, are the afterglow with full IGM absorption (dotted 
curve), and attempts to reproduce this profile with a damped Lya absorption system 
in the host galaxy (dashed curves). (Note, however, that damped absorption of this 
type could be suppressed by the ionizing effect of the afterglow UV radiation on 
the surrounding interstellar medium of its host galaxy|5Hj.) Most importantly, the 
overall spectral shape of the Lya trough carries precious information about the 
neutral fraction of the IGM at the source redshift; averaging over an ensemble of 
sources with similar redshifts can reduce ambiguities in the interpretation of each 
case due to particular local effects. 
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Fig. 3. Redshift distribution of all GRBs as compared to that measured by flux- 
limited surveys (taken from [5]). (a) Fraction of bursts that originate at a redshift 
higher than z vs. z. The data points reflect ~ 20 observed redshifts. (6) Fraction 
of bursts per logarithmic interval of (1 + z) vs. z. Solid lines: All GRBs for star 
formation through atomic line cooling. Dot-dashed lines: Expected distribution for 
Swift. Long-dashed lines: Expected distribution for BATSE. Note that the curves 
for the two flux-limited surveys are very uncertain due to the poorly-determined 
GRB luminosity function. Nevertheless, GRB number counts could provide a rough 
measure of the cosmic formation rate of massive stars at high redshifts. 



2.3 Statistics 



The upcoming Swift satellite (see http://swift.gsfc.nasa.gov/), planned for 
launch by the end of 2003, is expected to detect about a hundred GRBs per 
year. Figure 3 shows the expected redshift distribution of all GRBs (solid 
lines), under the assumption that the GRB rate is simply proportional to the 
star formation rate UJ. The long-dashed and the dot-dashed lines show the 
same rate, but taking into account the limited sensitivity of the BATSE and 
Swift detectors, respectively. This calculation implies that about a quarter of 
all GRBs detected by Swift should originate at z > 5. This estimate is rather 
uncertain because of the poorly determined GRB luminosity function 
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Fig. 4. Predicted number of line-of-sight intersections of metal-enriched bubbles 
around the first galaxies per redshift interval at a redshift z — 8 (taken from |16p. 
This statistic was calculated as a function of the equivalent width threshold of the 
metal absorption lines produced by these bubbles, for typical parameter choices of 
the theoretical model (see |16| for details). Thick solid, long-dashed, short-dashed, 
and dotted curves are for C II, O I, Si II, and Fe II metal lines, respectively. Thin 
solid and short-dashed curves are for C IV and Si IV, respectively. 

In principle, the rate of high-redshift GRBs may be significantly suppressed 
if the early massive stars fail to launch a relativistic outflow. This is possible, 
since metal-free stars may experience negligible mass loss before exploding 
as a supernova. They would then retain their massive hydrogen envelope, 
and any relativistic jet might be quenched before escaping the star. However, 
localized metal enrichment is expected to occur rapidly (on a timescale much 
shorter than the age of the then-young universe) due to starbursts in the first 
galaxies and so even the second generation of star formation could occur in an 
interstellar medium with a significant metal content, resulting in massive stars 
that resemble more closely the counterparts of low-redshift GRB progenitors. 

We note that the known population of optically-dark GRBs [22| is most 
likely associated with dust-obscured GRBs at low redshifts, although a small 
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fraction of these bursts could be optically-dark due to Lya absorption by the 
IGM at high-redshift. 

2.4 Progenitors 

In the previous section, we have assumed for simplicity a constant efficiency 
of forming GRBs per unit mass of stars. This simplifying assumption may 
either overestimate or understimate the frequency of GRBs. Metal- free stars 
are thought to be massive ^ |S] and their extended envelopes may suppress 
central jets within them (which may be produced through the collapse of their 
core to a spinning black hole). On the other hand, low-metallicity stars are 
expected to have weak winds with little angular momentum loss during their 
evolution, and so they may preferentially yield rotating central configurations 
that produce GRB jets after core collapse |21U17| . The basic question: should 
GRBs be common in the first generation of metal-free stars ? remains open for 
further study. 

3 Swift: A Challenge for Infrared Follow-Ups 

A small fraction (~ 10-25%) of the GRB afterglows detected by Swift are 
expected to originate at redshifts z > 5. This subset of afterglows can be 
selected photometrically using a small telescope, based on the Lya break at a 
wavelength of 1.216 /mi [(l + z)/10], caused by intergalactic HI absorption. The 
challenge in the upcoming years would be to follow-up on these candidates 
spectroscopically, using a large (10-meter class) telescope. A high-resolution 
spectrum can then be used to trace the ionization state (see Fig. 2) and metal 
enrichment state of the gas along the line-of-sight. For example, Figure 4 
illustrates the expected number of intersections of enriched bubbles of super- 
nova ejecta from early galaxies, as a function of the equivalent width of the 
metal absorption lines they produce (see JHj for details) . Measurement of this 
statistic for different metals can be used to constrain the initial mass function 
(through element abundances) and the mechanical energy output (through 
bubble sizes) of the first stars, as well as the significance of molecular H2 cool- 
ing in forming stars within the first mini-halos (through the number density 
of small bubbles). 

Based on the arguments mentioned in §1, GRB afterglows are likely to 
revolutionize observational cosmology and replace traditional sources like 
quasars, as probes of the IGM at z > 5. We caution, however, that this is 
a theoretical expectation (the highest redshift GRB detected so far j2j is at 
z ~ 4.5) and as such it is subject to uncertainties about the production of 
relativistic outflows by the first generation of stars. 

However, note that even if GRB outflows are quenched in the massive 
envelopes of low-metallicity stars, one would expect them to appear at high 
redshifts. This follows from the fact that metallicity obtains a large range 
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of possible values (in different spatial regions) at any slice of cosmic time. 
For example, the broad emission lines of all known quasars show evidence 
for high metallicity gas |2(J1 I14| . and so it is clear that there are pockets 
of highly enriched gas near the centers of massive galaxies at z > 6, where 
star formation reached an advanced stage and GRB production should be as 
likely as it is in the local universe. The cores of the most massive galaxies at 
z > 6 simply predate the characteristic evolution of the rest of the universe 
at a much lower redshift, because these cores reside in unusually overdense 
cosmological environments |24| . Such environments are likely to produce GRB 
progenitors in the same way that less overdense cosmological regions produce 
them at a lower redshift. Swift has the potential to educate us whether GRBs 
are produced only within these massive galaxies or in the more common sites 
of star formation at early cosmic times, namely the low-metallicity dwarf 
galaxies. 

In either case, the near future promises to be exciting for GRB astronomy 
as well as for studies of the high-redshift universe. 
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